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ABSTRACT 

Recent research has pointed to the ubiquity and abundance of between-generation epigenetic inheri¬ 
tance. This research has implications for assessing disease risk and the responses to ecological stresses and 
also for understanding evolutionary dynamics. An important step toward a general evaluation of these 
implications is the identification and estimation of the amount of heritable, epigenetic variation in popu¬ 
lations. While methods for modeling the phenotypic heritable variance contributed by culture have 
already been developed, there are no comparable methods for nonbehavioral epigenetic inheritance 
systems. By introducing a model that takes epigenetic transmissibility (the probability of transmission of 
ancestral phenotypes) and environmental induction into account, we provide novel expressions for 
covariances between relatives. We have combined a classical quantitative genetics approach with informa¬ 
tion about the number of opportunities for epigenetic reset between generations and assumptions about 
environmental induction to estimate the heritable epigenetic variance and epigenetic transmissibility for 
both asexual and sexual populations. This assists us in the identification of phenotypes and populations 
in which epigenetic transmission occurs and enables a preliminary quantification of their transmissibility, 
which could then be followed by genomewide association and QTL studies. 


E PIGENETIC inheritance involves the transgenera- 
tional transmission of phenotypic variation by means 
other than the transmission of DNAsequence variations. 
Cellular epigenetic inheritance, where transmission of 
phenotypic variation involves passing through a single¬ 
cell stage (the gametic stage in sexually reproducing 
multicellualr organisms), is now recognized to be an im¬ 
portant and ubiquitous phenomenon and the mecha¬ 
nisms underlying it are becoming elucidated (Jablonka 
and Lamb 2005; Allis et al. 2007; Jablonka and Raz 
2009). Epigenetic inheritance occurs between gener¬ 
ations of asexually and sexually reproducing organisms, 
directly affecting the hereditary structure of populations 
and providing a potential mechanism for their evolution 
(Jablonka and Lamb 1995, 2005; Bonduriansky and 
Day 2009; Jablonka and Raz 2009; Verhoeven et al. 
2009). It is therefore necessary to develop tools to study 
its prevalence and estimate its contribution to the heri¬ 
table variance in the population (Bossdore et al. 2008; 
Johannes et al. 2008, 2009; Richards 2008; Reinders 
et al. 2009; Teixeira et al. 2009). 

Unlike epigenetic inheritance, the inheritance of 
cultural practices in human populations has received a 
great deal of theoretical attention. Models of cultural 
inheritance and of interacting cultural and genetic effects 
have been suggested (Cavalli-Sforza and Feldman 
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1973; Rao et al. 1976; Cloninger et al. 1978; Boyd and 
Richerson 1985; Richerson and Boyd 2005). These 
models study the effects of cultural transmission and 
analyze the way in which it affects the distribution of 
cultural practices in the population. Other aspects of 
transgenerational effects are revealed through the study 
of maternal or indirect genetic effects (Kirkpatrick 
and Lande 1989; Wolf et al. 1998) and transgenera¬ 
tional genetic and epistatic effects within the context of 
the “missing heritability” problem (Nadeau 2009). 

The simple models described in this article focus on 
the transmissibility of epigenetic variations rather than 
on the magnitude of the phenotypic expression. In that 
respect, epigenetic inheritance is generally simpler to 
model than cultural inheritance since it commonly 
involves only vertical transmission (from parent to off¬ 
spring) . Crucially, during early development, as well as 
during gametogenesis and meiosis, some of the parental 
epigenetic information is restructured and reset. It is 
therefore necessary to explicitly include in models of 
epigenetic inheritance the number of developmental- 
reset generations. The number of developmental-reset 
generations between relatives may differ even when 
genetic relatedness is the same: for example, the re¬ 
latedness between parent and offspring is 0.5 and so 
is the relatedness between sibs (on average), but the 
number of developmental-reset generations is one and 
two, respectively. These considerations are also valid for 
asexual organisms, if it is assumed that some form of 
reset occurs during the cell cycle between divisions. In 
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this case we can test the models and measure the con¬ 
tribution of epigenetic inheritance in well-defined ex¬ 
perimental conditions, in pure lines. 

To estimate the amount of heritable epigenetic 
variation, we need to define several concepts: heritable 
epigenetic variability, the reset coefficient, and its com¬ 
plement, the epigenetic transmission coefficient. Heri¬ 
table epigenetic variability refers to phenotypic variability 
that is determined by epigenetic states that are environ¬ 
mentally induced and also possibly inherited from pre¬ 
vious generations. The heritable variations on which 
epigenetic heritability depends can arise spontaneously 
(as a particular type of developmental noise), or they 
can be environmentally induced. Once present, these 
variations can be vertically transmitted. For example, 
variations in methylation patterns between individuals 
may contribute to phenotypic variability even if these 
individuals are all genotypically identical. Such varia¬ 
tions have been found in several systems (Jablonka and 
Raz 2009). When the methylation marks are trans¬ 
mitted between generations, this will contribute to 
inherited epigenetic variability. The reset coefficient (v) 
refers to the probability of changing the epigenetic 
state during gametogenesis and/or early development, 
so that the new generation can respond to the present 
environmental conditions with no memory of past en¬ 
vironments. The complement of this reset coefficient, 
1 — v, is the coefficient of epigenetic transmissibility : the 
probability of transmitting the epigenetic state to the 
next generation, without reset. The epigenetic trans¬ 
mission coefficient should be taken as an abstraction 
that encapsulates all the potentialities of epigenetic 
inheritance related to the target phenotype, over a 
single generation. While we express the coefficient as a 
probability, it may also be interpreted as ratio coefficient, 
representing the portion of the epigenetic value that 
is transmitted to the next generation. In terms of the 
model these two interpretations are equivalent. The 
model also makes the simplifying assumption that the 
reset coefficient is a constant of the population, al¬ 
though in reality it might assume some distribution. 

We determine the above terms for both asexually and 
sexually reproducing organisms within a quantitative 
genetics framework. We suggest that terms including 
the effects of reset and epigenetic transmissibility be 
included in the classical model of phenotypic variance, 
showing how these terms can be calculated from the 
covariances between relatives. Morphological and phys¬ 
iological phenotypes can be assessed in this way, and 
if the results point to an epigenetic component of in¬ 
heritance, the underlying epigenomic bias can be 
investigated by employing QTL and association studies 
(Johannes et al. 2008; Reinders et al. 2009). A recent 
study using these methods provides evidence that epi¬ 
genetic variation can contribute significantly to the 
heritability of complex traits (~30% heritability), in¬ 
troducing transgenerational stability of epialleles into 


quantitative genetic analysis (Johannes etal. 2009). Our 
model, which entails direct measurements of pheno¬ 
typic covariation, should provide only preliminary and 
rough approximations of epigenetic transmissibility for 
target traits and populations. 

MODELS AND MEASURES OF EPIGENETIC 
VARIANCE, RESET, AND TRANSMISSIBILITY 

We consider a continuous trait with genetic, epige¬ 
netic, and environmental variability. First we develop 
the model for the case of asexual reproduction; sexual 
outbreeding populations follow. The phenotypic value 
of an individual ( P ) is given by the sum of the genetic 
value (G), the heritable epigenetic contribution (C), 
and the environmental deviation (£). Instead of the 
standard quantitative genetics model P = G + E we 
use P = G + C + E as our basic model, assuming 
independence and additivity, for simplicity. The two 
models are effectively equivalent, since in the classical 
equation E includes all nongenetic effects, of which 
our C may be one component. Alternatively, C may be 
included in G, since the variance attributable to herita¬ 
ble effects is sometimes assumed to be purely of genetic 
origin. The epigenetic value of an offspring depends 
on whether its epigenetic state was reset. If reset, the 
epigenetic contribution is determined by the inducing 
environment, I t . Offspring that are not reset inherit the 
epigenetic state, C, of their parents. An offspring resets 
its epigenetic state with probability v (the reset co¬ 
efficient) and inherits the parent’s epigenetic state 
with probability 1 — v (the epigenetic transmissibility 
coefficient). 

The inducing environment, I t , can be an exposure 
to an environmental signal or stress ( e.g ., a heat shock, 
an exposure to a chemical such as the demethylating 
agent 5-azacytidine, or exposure to an androgen sup¬ 
pressor such a vinclozolin), which elicits a developmen¬ 
tal reaction. Alternatively, it can be a mutation that has 
developmental effects that persist after the mutation 
itself has been segregated away. An example is the effect 
of ddml mutation in Arabidopsis, which leads to wide- 
ranging demethylation, with some demethylated pat¬ 
terns persisting many generations after the original 
mutation has been segregated away and the normal 
wild-type allele has been introduced (Reinders et al. 
2009; Teixeira et al. 2009). Another type of genomic 
shock leading to heritable epigenetic responses is hy¬ 
bridization followed by polyploidization (allopolyploidy). 
In these cases, notably in the case of Spartina hybrids, 
the newly formed allopolyploids acquire genomewide 
epigenetic variations, some of which are heritable for 
many generations (Salmon et al. 2005). When the in¬ 
ducing environment involves an environmental change, 
random or periodic changes in the inducing environ¬ 
ment maintain epigenetic variability in the population 
and can render a selective advantage to epigenetic 
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inheritance compared to noninducible genetic or 
nonheritable plastic strategies (Jablonka et al. 1995; 
Lachmann and Jablonka 1996; Ginsburg and 
Jablonka 2009). 

The epigenetic contribution of a reset offspring in 
year t is a random value C* taken from a distribution 
determined by the current inducing environment I ,; it is 
independent of the parental epigenetic state, C. Due 
to the temporal variation and epigenetic transmissibil- 
ity, the population mean of the epigenetic contribu¬ 
tion changes i n tim e as an autoregressive process, Y, = 
(1 — rjf, | + vC'i' 1 . In contrast to the inducing environ¬ 
ment, the microenvironmental factors that determine 
the environmental deviation ( E) vary among individuals 
within a generation, but on average do not change in 
time, such that the expected value of E is zero. We 
assume that G, C, and E are statistically independent 
such that 

Vp=V G +V c +V E . (1) 


during this process that developmental reset may occur 
(there may be a degree of reset in the parent too, in 
some cases). 

We begin by describing the specifics of the mathe¬ 
matical model as it pertains to asymmetric asexual 
reproduction, follow up with the unicellular symmetri¬ 
cal asexual case, and end with the sexual reproduction 
case. 

Asexual asymmetric reproduction: Consider a par¬ 
ent with phenotype P = G + C + E undergoing 
asexual reproduction. The mean phenotype of its 
offspring is 

Po= G+{l~v)C+vd r] . (2) 

Note that there is no microenvironmental term in (2), 
since we employed the mean offspring value, and the 
mean microenvironmental deviation is defined to be 
zero across the population. 

The parent-offspring covariance 


Two questions arise when the epigenetic term is in¬ 
troduced: (1) Can we distinguish, by quantitative genetic 
analysis, between heritable epigenetic and genetic in¬ 
heritance? and (2) How can we estimate the epigenetic 
transmissibility coefficient, (1 — v)? 

Both questions can be addressed by comparing co- 
variances between relatives. There are three unknowns 
to be estimated, the genetic variance, Vq, the heritable 
epigenetic variance V G , and the epigenetic transmissi- 
bility coefficient (1 — v ), since these are the only param¬ 
eters that contribute to the covariances. Therefore we 
need to compare three covariances between relatives 
that contain the genetic and epigenetic variance 
components in different proportions. The simplest 
way is to measure the covariance between parents and 
offspring, between sibs, and between uncles and neph¬ 
ews. We employ these anthropomorphic relational 
terms to maintain standard parlance within quantitative 
genetics. 

We differentiate between two major types of asexual 
reproduction: (i) asymmetric transmission, where a par¬ 
ent organism is not identical to its offspring following 
cell division (this is the case following reproduction 
by budding, with the parent organism able to bud off 
several offspring over time), and (ii) symmetrical cell 
division, where the parent cell is also the offspring cell, 
as occurs in many unicellular organisms and in cell 
lineages. In epigenetic terms, symmetrical reproduction 
leads to common reset potential for both daughter cells 
(either both reset or neither of them do), but the 
inducing random epigenetic contributions they obtain 
are independent. In this case, single cells have to be 
followed over generations for the family relations to 
be established. We assume that asexual asymmetric re¬ 
production involves a developmental process of the 
offspring that the parent need not go through, and it is 


COVop = COV[(G+ C + E), (G + (1 
= V G + (1 -v)V c 


W )C+wC f (r) )] 


( 3 ) 


follows directly from the assumption that G, E, C, and 
Ct'' 1 are mutually independent and from the fact that the 
covariance of any individual with the mean value of a 
number of relatives is equal to its covariance with any 
one of those relatives (Falconer and Mackay 1996). 
The covariance of sibs equals the variance of the true 
means of sib groups, which is 


COVsib 


Var 


G+ (1 


v)C+vCi r) 


= Vb + (1 -vfVc, 


( 4 ) 


assuming C\ ’ : does not covary between siblings. 

Note that the heritable epigenetic variance contrib¬ 
utes less to the covariance between sibs than to the 
covariance between parents and offspring. Under asex¬ 
ual reproduction, purely genetic variability results in 
equal covariances between relatives. A significant differ¬ 
ence between COVop and COVsib therefore may sug¬ 
gest heritable epigenetic variability. 

The covariance between uncles and nephews is (see 
appendix) 


COV un . = V G + (1 - vfV c . (5) 

Figure 1, A and B, depicts how the epigenetic con¬ 
tribution to the similarities between relatives depends 
on the opportunities for epigenetic resets, in asexual 
species, for parent-offspring, between-sibling, and uncle- 
nephew relations. 

From Equations 3-5 one can calculate the variance 
components and the epigenetic transmissibility coeffi¬ 
cient as 
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1 _ COVsib — COVun 

COVop — COVsib 
COVop — COVsib 
C ~ v(l - v) 

V G = COVop - (1 - v) V c . (6) 

If the covariance estimates COVsib and COVop are 
approximately equal, then it can be shown from the 
derivation of (6) that either V (: is near zero (i.e., a single 
value to the epigenetic effect) or the reset coefficient v 
is close to its extreme values of one or zero. We would 
not be able to distinguish between these three distinct 
possibilities, all of which are theoretically viable. This 
introduces a singularity point inherent in our model, 
which we acknowledge, and in such a case the model 
is inapplicable for the target phenotype and popula¬ 
tion. The same cautionary note applies to (8) below for 
COVsib = COV UN and to (9) for COV OP = 2 COV HS - 
Furthermore, v = 0 would mean that epigenetic effects 
become indistinguishable from genetic effects; simi¬ 
larly, v = I means that epigenetic effects are always reset 
and therefore indistinguishable from environmental 
effects. 

Note that we could extract estimates of heritability 
from our model, since it also generates Vq. Checking 
such estimates with reference to various traits against 
known heritability values should provide a further 
indication as to the soundness of our model. More¬ 
over, V c wi 11 be effectively zero if the population has a 
single recent ancestor (not withstanding the possible 
mitotic mutations over generations). In that case we 
should expect COVop to be equal to (1 — v) Vc, and any 
deviation from this equality may indicate a flaw in the 
model. 

As a consequence of the changes in the inducing 
environment, I t , across generations and consequently 
in the distribution of C\ r \ the epigenetic variance 
may change in time. Therefore the ancestors of the 
measured relatives, i.e., the parent in Equation 3, the 
parent of the sibs in Equation 4, and the parent of 
the uncle (grandparent of the nephew) in Equation 5, 
must belong to the same generation, and Vq is the 
variance of the epigenetic contribution in this ancestral 
generation. 

Asexual symmetric reproduction: Here the focus is 
on a population of unicellular organisms or cell line¬ 
ages. This case is naturally limited in lending itself to 
large-scale statistical quantitative analysis, as meticulous 
tracing of cell lineages across generations is required for 
producing the covariances. Since the offspring resets are 
not independent, we require examining also the co- 
variance of cousins, labeled COVcous> when solving 
for the three unknowns. Additionally, the possibility of 
a fluctuating inducing environment across generations 
introduces an extra term in three of the covariances (see 
appendix) , such that 


(single) Parent* 


Genetic relatedness = 1 

One chance for reset: 
epigenetic contribution 
proportional to (1-v) 



Offspring 1 
*asymmetric cell division 


Two independent 
chances for reset: 
epigenetic contribution 
proportional to (1-v) 2 


Offspring 2 


B (single) Grandparent* 



v 


Uncle 


Parent* 


Genetic 
relatedness = 1 


Three independent 
chances for reset: 
epigenetic contribution 
proportional to (1-v) 3 


*asymmetric cell division 


Chance for 
reset, v 



Offspring 


Figure 1.— (A) Parent-offspring and full sibling similarities 
in asymmetric asexual reproduction: the genetic relatedness is 
the same in both cases, while the epigenetic contribution dif¬ 
fers. (B) Uncle-nephew similarity in asymmetric asexual repro¬ 
duction: the epigenetic contribution is proportional to (1 — v) 3 . 

COVop = V G + (l-v)V c 
COVsib = V G + (1 - v)V c + w(l - v)(cj r) - rVi) 2 
COV UN = V G + (1 - v) 2 V c + v(l - v) 2 (c[ r) - r,_i) 2 
COVcous = V G + (1 - vf V c + i/(l - v) 3 (C, (r) - IVi) 2 . 

(7) 

Figure 2, A and B, depicts how the epigenetic con¬ 
tribution to the similarities between relatives depends 
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on the opportunities for epigenetic resets across single 
and multiple generations of offspring cells in an asexual 
species with symmetric reproduction. 

We then arrive at the solution for our three unknowns, 


_ COVun — COVcous 

COVsib — cov UN 

u COVop + (1 — w)COVsib — COVun 
V G = COVop - (1 - v) V c . (8) 


Sexual reproduction: We assume that the parents con¬ 
tribute equally to the epigenetic contribution of their 
nonreset offspring. This assumption is valid for chro- 
mosomally transmitted epigenetic information, such as 
methylation patterns, histone modifications, nonhis¬ 
tone DNA binding proteins, and nuclear RNAs. When 
heritable epigenetic variations are cytoplasmic or corti¬ 
cal, the mother alone transmits the epigenetic variation, 
and when only the female lineage is considered, the 
variations produced by these epigenetic inheritance 
systems can be analyzed in the same way as variations 
transmitted by asexual, asymmetrically reproducing 
organisms. Here we consider chromosomally transmis¬ 
sible epigenetic variations as well as variations mediated 
through the RNAi system and assume equal contribu¬ 
tion by both parents. It is likely that during sexual 
reproduction, the reset mechanism is more compre¬ 
hensive than in asexual reproduction. However, as in the 
asymmetric asexual scenario, the offspring here reset 
independently. 

The heritable epigenetic component creates several 
discrepancies from the standard quantitative genetic re¬ 
lations. Due to the larger fraction of epigenetic variance 
contributed to the parent-offspring covariance, COVop 
is more than twice the covariance between half sibs 
(COV 0 p > 2 COV HS ); however, this difference could 
also be explained by the effect of additive-additive 
interaction. A more conclusive comparison is between 
the uncle-nephew and the half-sib covariance. Uncles 
and nephews as well as half sibs are second-degree 
relatives with a single common ancestor, and therefore 
their covariance contains exactly the same genetic vari¬ 
ance components. However, the epigenetic variance 
contributes more to the covariance between half sibs 
(COV H s > COV UN ). Detecting a greater covariance 
between half sibs therefore would provide good evi¬ 
dence for epigenetic inheritance. 

Comparison of the full-sib and parent-offspring 
covariances might also reveal the presence of heritable 
epigenetic variation. The epigenetic contribution in¬ 
creases COVop more than COV F $ while the additive 
genetic contribution is the same; thus COV FS < COV Q p 
would be expected. However, dominance and maternal 
effects can increase COV F s considerably and hence 


Parent Cell* 



Single chance for 
reset, v, at parent cell 


Generation AM 
Generation N 

= 1 X 

Two dependent ' ' - 
chances for reset: 
epigenetic contribution 

Offspring Cell proportional to (1-v) Offspring Cell 
*symmetric cell division 

B Grandparent Cell* 


Single chance 
for reset, v 


Genetic relatedness = 1 
for all cell divisions 


Generation A 1-1 


Generation N 



Uncle Cell* 


Cousin Cell 
*symmetric cell division 


Generation 

N+1 


Between Cousins: I 
Three independent ' 
chances 

for reset: (1-v) 3 Nephew/Cousin Cell 


Figure 2.—(A) Parent-offspring and between-sibling sim¬ 
ilarities in symmetrical asexual reproduction: single chance 
for reset for both offspring. (B) Uncle-nephew and between- 
cousin similarities in symmetrical asexual reproduction: the 
parent cell and the uncle cell share a single chance for reset. 


somewhat mask the effect of heritable epigenetic vari¬ 
ation. The presence and magnitude of the obscuring 
dominance and maternal effects can be assessed by sib 
analysis (Falconer and Mackay 1996, p. 166). Sib 
analysis makes use of the fact that without dominance 
and maternal effects the full-sib covariance should equal 
twice the half-sib covariance. This remains true even 
in the presence of heritable epigenetic variation, so sib 
analysis can be applied to measure dominance and 
maternal effects without any distortion due to the her¬ 
itable epigenetic effects. Removing the dominance and 
maternal effects from the full-sib covariance facilitates 
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the detection of heritable epigenetic variability by 
comparison of the full-sib covariance with the parent¬ 
offspring covariance. For simplicity, we assume no domi¬ 
nance or epistatic effects [for instance, one-eighth of the 
additive-additive variance contributes to COV 0 p — 2 
COV H s in (9)] and suggest using paternal half siblings 
where maternal effects are absent. 

Figure 3, A and B, depicts how the epigenetic con¬ 
tribution to the similarities between relatives depends 
on the opportunities for epigenetic resets in sexual 
species. 

The variance components and the epigenetic trans- 
missibility coefficient can be estimated from the mea¬ 
sured covariances as 


Parent* 


Genetic relatedness = !4 

One chance for reset: 
epigenetic contribution 
proportional to (1-v) 


Chance for reset, v 



Offspring 1 proportional to (1-v) Offspring 2 


1 _ _ 2(COV H s — COVun) 

COVop — 2 • COVhs 
_ 2(COV OP - 2 • COVhs) 

C w(l - v) 

Va = 2-COV op —(1 -v)V c (9) 

(see appendix), where the uncle and the parent of the 
nephew are full sibs, and, as before, the parent in 
COVop, the parent of the half sibs, and the parent of the 
uncle (grandparent of the nephew) belong to the same 
generation. 

DISCUSSION AND CONCLUSION 

Transgenerational, ecologically or developmentally 
induced phenotypic variations have been studied 
mainly in the context of maternal effects and have been 
treated as temporally extended developmental effects 
(Wolf et al. 1998). The study of cultural transmission, 
on the other hand, was confined to studying the com¬ 
plex system of transmission of cultural practices (Boyd 
and Richerson 1985; Feldman and Laland 1996). As a 
consequence, quantitative studies of the transmission 
rather than the expression of heritable epigenetic vari¬ 
ations have only in recent years begun to surface, 
although molecular studies have been around for some 
time showing that epigenetic marks such as DNA meth- 
ylation patterns, protein marks, and RNA-mediated 
silencing can be inherited (reviewed in Jablonka and 
Raz 2009). 

In the absence of detailed molecular studies, it is 
often very difficult to detect epigenetic inheritance, 
since deviations from classical Mendelian ratios can 
be always explained by assuming interactions with 
modifier genes. In asexually reproducing, multicellular 
organisms, rapid and heritable phenotypic switches are 
usually explained within the framework of somatic 
mutations and somatic selection or various types of 
phase variations. 

It is instructive to suggest a formulation of the extent 
of heritable epigenetic variation. The epigenetic heritabil- 
ity, denoted here by y 2 , is the proportion of the total 


"Sexual: Showing only one of the two parents 


B Grandparent* 


v 


Uncle 


Parent* 



Genetic 

relatedness = % 


Three independent 
chances for reset: 
epigenetic contribution 
proportional to (1-v) 3 


* Sexual: Showing only one of 
the two parents, in each case 


Chance for 
reset, v 



Offspring 


Figure 3.— (A) Parent-offspring and full sibling similari¬ 
ties in sexual reproduction: two independent chances for re¬ 
set between offspring, while their genetic relatedness is g. (B) 
Uncle-nephew similarity in sexual reproduction: three inde¬ 
pendent chances for reset between uncle and nephew, while 
the genetic relatedness is |. 

phenotypic variance V P attributable to the potentially 
heritable epigenetic variance Vc, 


9 

r = 


Vc 

v P ‘ 


( 10 ) 


Contrary to the standard notion of heritability where 
V G is always heritable variance, the epigenetic herita- 
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bility, by definition, does not directly depend on the 
epigenetic transmissibility coefficient, (1 — w), but only 
on Vc■ However, the extent to which epigenetic herita- 
bility contributes to the measurable regression or cor¬ 
relation coefficients does depend on the epigenetic 
transmissibility of the epigenetic inheritance system. 
For example, in parent-offspring regression, a fraction 
(1 — v) of the epigenetic heritability always increases the 
regression slope: 


b Q p — 


1 Va+(1-p)Vc 

2 V P 


i[/t 2 + (i - w )r]. 


(ii) 


In general, the epigenetic contribution to the re¬ 
gression of two related individuals would be (1 — v) n y 2 , 
where nis the number of independent opportunities for 
epigenetic reset in the lines of descent connecting the 
relatives. 

The inducing environment, I t , may determine the 
epigenetic contribution after reset in a deterministic 
fashion, such as under extreme stress conditions, in 
which case the distribution of C\ r) becomes a delta-peak. 
In a situation where I t ceases to vary at generation k, 
the reset individuals always acquire a new epigenetic 
contribution from the same distribution C lr> , so the 
distribution of C will converge on the distribution of 
C (r) and r n = C (r ) + (1 — — C (r )). In the ex¬ 

treme case when I, both is time independent and 
unequivocally determines the induced epigenetic con¬ 
tribution, the distribution of C converges to a delta- 
peak; no heritable epigenetic variation is present. 

In the model presented in this article we assumed 
that v, the reset coefficient, is independent of the 
number of generations that the lineage spent in a 
particular environment prior to the switch to a new 
inducing environment. However, it is possible that the 
number of generations a lineage spends in a given 
environment may alter the probability of a phenotypic 
switch. Transfer experiments, in which the number of 
generations spent in each environment is varied, and 
v and 1 — v are measured for each case, may uncover 
such phenomena. Other factors that may affect the 
stability of epigenetic inheritance are the number of 
generations a chromosome is transmitted through 
only one sex ( e.g ., a male) and the number of gen¬ 
erations that chromosomes are continuously trans¬ 
mitted though an aged parent (Jablonka and Lamb 
2005). 

What values of the epigenetic transmissibility coeffi¬ 
cient and epigenetic heritability should alert us to the 
possibility that part of the observed variance is due to 
epigenetic variation? Ideally, we want to have threshold 
values for each phenotype under consideration, so that 
only values obtained through the model that lie above it 
would provide sufficient confidence for further molec¬ 
ular inspection (see the appendix for a suggestion for 
employing a workable probabilistic threshold). In terms 


of the epigenetic transmissibility coefficient, it is reason¬ 
able to treat any value above zero as indicative of pos¬ 
sible heritable effects, even if only for a very small 
number of generations and individuals. 

The detection and measurement of epigenetic heri¬ 
tability may suggest new testable interpretations of a 
number of observations. For example, experiments mea¬ 
suring heritability in inbred lines showed that heritabil¬ 
ity increases fairly rapidly in inbred sublines of different 
organisms (Grewal 1962; Hoi-Sen 1972; Lande 1976). 
These results were interpreted to mean that many genes 
contribute to a character and that the rate of mutation 
of these genes is high. The mutation rates calculated 
from the data with the assumption that mutations alone 
contribute to heritability were two to three orders of 
magnitude higher than the traditionally estimated rate 
of classical mutation (Hoi-Sen 1972). But the high 
heritability appearing after relatively few generations of 
subline divergence can be reinterpreted as being partly 
due to heritable epigenetic variability, especially in 
view of more recent molecular studies that show that 
epigenetic inheritance occurs in pure lines (Johannes 
et al. 2008, 2009; Reinders et al. 2009). When a subline 
of a highly inbred line is started from a pair of indi¬ 
viduals, epigenetic variations (epimutations) as well as 
mutations will accumulate as the subline diverges. As 
the rate of epimutation by epigenetic reset is often 
much higher than that of mutation, heritable epige¬ 
netic variability will grow more rapidly than genetic 
variability. The fraction of epigenetic variability should 
be relatively high, especially during the first generations 
of a subline divergence, when it is the sole contributor to 
heritable variance. Since the epigenetic heritability can 
be measured from the covariances between relatives 
within the line, this hypothesis can be tested. 

The hypothesis that there is a heritable epigenetic 
component of phenotypic variance has the implication 
that selection can be effective in genetically pure lines. 
Some observations indicate that this is indeed the case 
(Brun 1965; Ruvinsky et al. 1983a,b; see Johannes 
et al. 2009 for a precise estimate based on QTL analysis). 
The response to selection in one generation is deter¬ 
mined by the midparent-offspring regression coeffi¬ 
cient, b = h~ + (1 — v)y 2 . Since pure lines may acquire 
epigenetic variability relatively quickly, and this varia¬ 
bility contributes to the regression, the response to 
selection enabled by the epigenetic heritability can be 
readily measurable. The realized heritability in a genet¬ 
ically pure line will then estimate (1 — v) y 2 . Note that 
once selection is relaxed, all effects achieved by previous 
selection on the epigenetic variability are sooner or later 
lost (at a rate depending on epigenetic transmissibility), 
unless the phenotype is by then stabilized by genetic 
assimilation (Ruden et al. 2005, 2008; Shorter and 
Lindquist 2005). 

We have shown that epigenetic heritability can be 
measured for both asexually and sexually reproducing 
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organisms. A promising system in which to measure 
epigenetic heritability is in inbred lines of plants, since 
plants seem to have particularly sophisticated cellular 
mechanisms of adaptation, and several cases of epige¬ 
netic inheritance based on chromatin marks have been 
reported in plants (Jablonka and Raz 2009). The 
phenotypic effects of altering environmental conditions 
have been studied in plants, and data on the frequency of 
phenotypic change are known for some cases (Bossdorf 
et al. 2008; Johannes et al. 2008, 2009; Richards 2008; 
Reinders et al 2009; Whittle et al. 2009). 

The novelty of our approach lies in the attempt to 
quantify a measure of epigenetic variability using classi¬ 
cal notions of familial phenotypic covariances. Other 
approaches have looked at the actual epigenetic vari¬ 
ability, with the goal of producing a quantifiable mea¬ 
sure of the variability with respect to a certain locus or a 
related group of loci. One such research looked at 
epigenetic variability in the germline, formulating an 
average methylation-intensity vector for each locus/ 
individual, by looking at the sum of the methylated 
cytosines for each different cytosine position. The 
degree of epigenetic dissimilarity was then defined via 
the Euclidean distance between the vectors of the two 
individuals (Flanagan et al. 2006). 

A complementary approach to that presented here 
is proposed in Slatkin (2009). The model assumes that 
disease risk is affected by various diallelic genetic loci 
and epigenetic sites and takes allele frequencies and 
gain and loss rates of inherited epigenetic marks as input 
parameters. The mathematical iterative process of a two- 
state Markov chain is then utilized to simulate the state 
of the contributing loci after a number of generations. 
A central perspective that we believe is common to 
Slatkin’s approach and ours is precisely this introduc¬ 
tion of parameters for the quantitative effects of an 
inducing environment and reset coefficient and the 
realization that genetic and epigenetic factors contrib¬ 
ute differently to familial similarities (recurrence risk, in 
Slatkin’s case). However, our approach is distinct in that 
we employ and extend the standard quantitative genetic 
approach and utilize familial correlations in observed 
phenotypes as inputs to our model; it does not require 
prior estimations of gain and loss rates of inherited 
epigenetic marks or their frequency. 

Studies focusing on the timing of changing the en¬ 
vironment and applying an inducer can shed light on the 
sensitive periods of development, during which changed 
conditions can have long-term heritable effects. It has 
been shown that developmental timing is crucial for 
changed conditions to have long-term effects: for exam¬ 
ple, administering an androgen suppressor to a pregnant 
rat leads to heritable epigenetic effects only if the drug is 
administered during a sensitive period between days 8 
and 15 of pregnancy (Anway et al. 2005). 

Environmental insults, such as pollution, cause hu¬ 
man (and nonhuman) diseases that may carry over for a 


number of generations (Jablonka 2004; Gilbert and 
Epel 2009). For instance, early paternal smoking was 
associated with greater body mass index at 9 years of age 
in sons (Pembrey et al. 2006). Such carry-over effects 
may be uncovered by careful search for transgenera- 
tional effects. The ability to detect and quantify a 
preliminary rough measure of transgenerational, envi¬ 
ronmentally induced diseases could become part of 
epidemiological studies. 

The model presented here provides an initial esti¬ 
mate of epigenetic transmissibility that can be the basis 
for further molecular studies. We are well aware that a 
complete analysis of the developmental, and possibly 
evolutionary, effects of epigenetic inheritance has to 
include both transmissibility and expressivity and that 
QTL methodology is indispensable. However, a pre¬ 
liminary ability to empirically detect and quantify a 
rough measure of transgenerational epigenetic effects 
can give epidemiologically important information and 
assist in narrowing and directing the search domain for 
molecular epigenetic sequencing. 
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APPENDIX: DERIVATION OF FAMILIAL COVARIANCES 

Asexual asymmetric reproduction: For derivation of COV UN consider an asexual parent (P) with two offspring ( X 
and U) and an offspring of Xdenoted by N. U and N are uncle and nephew. The uncle was born in generation t— 1 and 
the nephew in generation t. Given the genotypic value, Gp, and the epigenetic contribution, Cp, of the common 
ancestor, the expected phenotypic value of the uncle is 

Pu= G p + (1-v)C p + vc¥ 1 . (Al) 

The nephew inherits the epigenetic contribution C P only if neither it nor its parent resets the epigenetic state. If the 
nephew is not reset but its parent is reset, then it inherits C t \; if the nephew resets, then it takes a random value Cj r \ 
Hence the expected phenotypic value of the nephew is given by 

Pn = Gp + (1 - v) 2 C P + v(l - v)c\ r \ + vC[ r) . (A2) 

From Equations Al and A2 the uncle-nephew covariance is 

COV UN = COV[Pu, P N ] = V a> + (1 - v) 3 V Ce . 


(A3) 
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Asexual symmetric reproduction: This case is complicated by the fact that due to simultaneous reset, a fluctuating 
inducing environment affects some of the familial covariances. The simultaneous reset of offspring does not affect the 
covariance between the parent and one of its offspring; hence 


COVqp = V G + (1 - v)V c 


(A4) 


just as in case of asexual asymmetric reproduction (see Equation 3 of the main text). To obtain the covariance between 
sibs, we must use conditional expectations and covariances. Within one family (with the parent having phenotype P= 
G + C + E), the conditional offspring mean is 


Po 


G + C _ if nonreset 

G + cF if reset 


(A5) 


so that we get COVsibs|nonreset = Var(G+ C) = V c + V c , COVsibsireset = Var(G+ Cj' 1 ) = V c [since we assumed 
that C\' ) does not covary between siblings]. By the law of total covariance, the unconditional covariance equals the 
mean of conditional covariances, 


(1 - v)(V G + V c ) + vV G = V G + (1 - v)V c 


plus the covariance of conditional means across the pop ulati on. The population mean of the nonreset offspring 
is G + r t _i, whereas the mean of reset offspring is G + C* . Define X to be the conditional mean of the first 
offspring in a family and Y the conditional mean of the second offspring in a family. Their bivariate distribution is 

Prob(X = Y = G + T,_j) = 1 — v (both nonreset); 

Prob(X = Y = G + cl r> ) = v (both reset). 


The covariance of con ditio nal means across the populatio n is t herefore COV(X, Y) = E(XY) — E(X)E( Y) = ( 1 — u ) 
(G + IVi) 2 + v(G + cl r) f — ((1 — v)(G + T,_i) + v(G + c\ r) )f, which readily simplifies to w(l — «)(cF — 
T^j) 2 . Putting all terms together, the unconditional covariance between sibs is given by 


COVsib = V G + (1-v)V c + v(l - v)(cl r) - IVO 2 . (A6) 

Note that if the inducing environment becomes constant across generations, then the autoregressive process T( = 
(1 — u)r,_! + vC <r > converges to T = C (r> so that the last term of Equation A6 vanishes, ft can be shown that unlike the 
covariance, the correlation of sibs increases with a large fluctuation in the mean of the inducing environment, but 
decreases with an increase in the variance of the inducing environment. A fluctuating inducing environment affects 
COVsib only in the symmetric case, since in the asymmetric case the bivariate distribution of siblings is such that the 
covariance of conditional means is always zero : 

Prob(X = Y = G + r ( _i) = (1 - v) 2 ; Prob(X =Y=G+ C, W ) = v 2 

Prob(X = G + r,_! and Y = G + Q (r) ) = (1 - v)v; 

Prob(X =G+ Cl r) and Y = G + T,^) = w(l - v). 

The covariance of conditional means is therefore 

COV(X, T) 

= E{XY) -E(X)E(Y) 

= (1 - vf{G + IVO 2 + u 2 (G + cF) 2 + 2v(l - v)(G + T ( _0(G + cF) 

- ((1 - u) 2 (G + r ( _i) + u 2 (G + cF) + (1 - v)v(G + r ( _i) + v(l - v)(G + gF )) 2 = o 

and similarly it is zero for the uncle-nephew covariance (asexual asymmetric case). 

For the calculation of uncle-nephew and cousin covariances we express the conditional expectations of the uncle U 
and of the nephew N, 

_ f G> + Cp if U nonreset 

U \ Gp + cF if U reset 


and 
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if U nonreset, N nonreset 
if U nonreset, N reset 
if U reset, N reset 
if U reset, N nonreset, 

where Gp is the genetic value and C P the epigenetic contribution of the parent of the uncle. In accordance with the 
respective probabilities for reset, the mean of conditional covariances is 

COV(Gp + C P , Gp + C P )(1 - v) 2 + COV(Gp + C P , Gp + C$)v(l - v) 

+ COV(Gp + c f (r) , Gp + c[ r X)v 2 + COV(Gp + C t (r) , Gp + C, (r) )u(l - v) 

= (1 - u) 2 (Vgp + V Cp ) + 1/(1 - vjVc, + 1/(1 -v)V Gp + V 2 V 0p = + (1 - v) 2 V C? . 

To obtain the unconditional covariance, we need to add the covariance of conditional means. If X is the conditional 
mean of the uncle and F is that of the nephew, we have the bivariate probability, 

Prob(X = F = G + IV0 = (1 - v) 2 ; Prob(X = G + T,_i and Y=G+ ) = w(l - v) 

Prob(X =Y=G+ C, W ) = v{\ - v); Prob(X = G + C f (r) and Y = G + c\X ) = v 2 . 

The covariance of conditional means comes out to be u(l — w) 2 (C / (r) — iyi) 2 . For the unconditional covariance 
between uncle and nephew we get 

COV UN = V G + (1 - v) 2 V c + w(l - v) 2 {c[ r) - r ( _!) 2 . (A7) 

For the between-cousin covariance, we first calculate the mean of conditional covariances, 

COV(Gp + C P , Gp + C P )(1 - w) 3 + COV(Gp + C P , Gp + C$)t/(1 - v) 2 

+ COV(Gp + C\ r } x , Gp + C P )u(l - v) 2 + COV(Gp + C$, Gp + C$)(l - v)v 2 
+ COV(Gp + Q (r) , Gp + c\ r) )v{ 1 - v) 2 + COV(Gp + C t (r) , Gp + C ( ( ^)t/ 2 (1 - v) 

+ COV(Gp + c£\, Gp + C[ r) )v 2 { 1 - v) + COV(Gp + d+i, Gp + C^u 3 
= Top + (1 - v) S Vc P 

and then add the covariance of conditional means. We denote X as one cousin and Fas the other cousin (where U is 
the father of the first cousin) and express their bivariate probability: 

Prob(X =Y=G + T ( _i) = (1 - v) 3 ; Prob(X = G + F,_, and F = G + C$) = (1 - v) 2 v 
Prob(X = G + C\X and F = G + r t _i) = (1 - vfv; 

Prob(X =Y=G+ C ( t r> ) = i/(l - v) 2 ; Prob(X =G+ c[ r) and F = G + ) = w 2 (l - v) 

Prob(X = G+ c£\ and F = G + C f (r) ) = u 2 (l - v); Prob(X =Y = G+ C l fX ) = v 2 . 

The covariance of conditional means comes out to be u(l — v) 3 (Ct — r^) 2 . So we have 

COVcous = V G + (1 - vfV c + t/(l - n) 3 (d r) - T^) 2 . (A8) 

From (A4), (A6), (A7), and (A8) we solve for our three unknowns as expressed in (8), where the unknown term 
( cj r) — r,_!) 2 cancels out. 

Sexual reproduction: Here A denotes the breeding value so that Va is the additive genetic variance; the mean 
breeding value is scaled to zero. With paternal phenotype P= A + C + £the offspring mean phenotype under random 
mating for paternal half siblings is 

Po = ^ + ^(1 - v)C + ^(l - w)r t _! + vC\ r) . 

The offspring-parent and half sibling covariances are 


Pn = 


Gp + Cp 
Gp+Cgl 
Gp + C\ r \ 


Gp + c; 


■M 
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COVop = - Va + — (1 — v) Vc\ 

COVhs = Var(7b) = \ V A + \ (1 - vf V c . (A9) 


For the derivation of uncle-nephew covariance (the parent of the uncle should come from the same generation as 
the parents above), let M and F denote a mother and a father with two offspring, X and U, and an offspring of X 
denoted by N; U and N are uncle and nephew as before (the other parent of the nephew is an individual taken 
randomly from the population). The expected phenotypic values of the uncle and his sib X are the same. 


— — Am T A] 

Pi j = Px = „ 


+ (1 - v) 


Cm + Cf 


+ vC 


(r) 
i-i • 


(A10) 


If the nephew did not reset his epigenetic state, then he inherits half the epigenetic contribution of X plus half the 
epigenetic contribution of a randomly chosen individual in generation t — 1. Therefore the nephew’s expected 
phenotypic value is 


Pn = ^A x + * (1 - v) C x + * (1 - + vc\ r) 

(r) 

Am + Ap . >.2 Cm + Cf , -.Q-i 1 n , r „(r) 

= —^— + (!-«)—^— + w(i-w)— +-(i - ior ( _t + vc\ 


(All) 


From Equations A10 and A11 the uncle-nephew covariance is calculated as 

COV UN = COV(Pb, Ts) = \ V A + * (1 - vfV c , 


(A12) 


where V A = Var(A M ) =Var(A r ) and V (: = Var( G M ) = Var(C F ). 

Identifying classes of individuals with large epigenetic contributions: The objective is to find an expression for the 
probability that, for a given phenotypic deviation, P, \ C\ is greater than a certain threshold K, i.e., that this value of Cwill 
be at the extreme of its distribution, formally, that Prob (| C | > K \ P). The underlying assumption is that higher values 
of | C| generally indicate the presence of more epigenetic factors, separately identifiable. This perspective has the 
added benefit of directing the search to classes of individuals or pure lines corresponding to certain domains of 
phenotypic values, as the probability is conditional on P. 

We extend the model proposed in Tal (2009) to include our heritable epigenetic component, C, in compliance 
with P = G + C + E, and from statistical independence, Vp = V G + Vc+ V E . The only extra assumption we introduce 
here to comply with the model in the main text is that our quantitative trait, P, is approximately normally distributed. 
This allows us to infer the normality of G, C, and E and arrive at expressions for the joint conditional distribution of G 
and C and subsequently to the expression of probabilities. The three variables in our quantitative model represent 
deviations from their respective means (so that the mean of Pis also zero), and as a consequence we employ the 
absolute value of Gin the expression for the probability. Without loss of generality we assume that the variance of Pis 1 
(Pis standardized so that hr is the variance of Gand y 2 the variance of C). So we define 

P ~ N{ 0, 1), G - 1V(0, /i 2 ), C ~ N{ 0, y 2 ), E ~ N(0, 1 - Ir - y 2 ) 0</i 2 + y 2 < 1. (A13) 

We need to arrive at the joint conditional density function of Gand C given Pand then identify the portion of the 
distribution that satisfies |C| > K. We denote tp„(x) the probability density function of a normal random variable X 
with zero mean and variance v. Now Pinduces ajoint conditional distribution of Gand Cgiven /r and y 2 . Let us denote 
this by c). From first principles of conditional probability we have 

F p,h*,y*(g, c ) =Pc,c\p(g, c\p) 

= pG,C,p(g . G P) = PP\G,c(p\g, c) ' p G ,c{g ■ C ) 

Pp(P ) Pp(P) 

= PP\G,dP\g > C ) ■ PcXg) ■ Pc(c) 

Pp{P) 

= Pe{P -g-c)- pc(g ) • pc{c) 

pp(p) ' 1 


Note that the move from pp\G,c(p\g , c) to p E (p ~ g — c) is justified since P = G + C + E. Now, in terms ofq> v (x) we have 
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Heritability=0.7, and various phenotypic values, p 



Figure Al.—The graphs of M for Ir = 0.7. Depicted for various values of the phenotypic deviation (in standard deviations), y 2 
varies from 0 to 0.3. These are the probabilities that the epigenetic deviation of an individual |C| is at the extremes of its distri¬ 
bution, Prob(| C\ > K \ P), where K = 1 (in standard deviations of P). 


F p,h 2 ,y 2 (g’ c ) = 


<Pl ~ g ~ c) ' 9ft 2 (g) • <P-y* (c) 

9l 0 ) 


(A15) 


Explicitly substituting tp 7 ,(x) and simplifying terms, we get the bivariate normal form, 


F (g, c) = 


P = 


2TTCT1CT2X/(1 - p 2 ) 

— hy 

V(!-**)(!-7*) ’ 


(-l/(2'(l-p 2 ))) [((g- ( T,)/CTl) 2 -2p((g-p 1 )/CTlH(c-p.,)/!T2) + ((c-p. 2 )/o 2 ) 2 


(A16) 


cfi 


= h\J\ - tr 2 = 7V 1 - Yh hi = h 2 p, p 2 = y 2 p- 


We may now employ F to find our conditional probability, Prob (| C \ > K \ P ). This probability is in fact the integral 
of the distribution of Fin the domain that satisfies |C| > K, denoted by M pjp(y~). We therefore have 


Heritability=0.3, and various phenotypic values, p 



Figure A2.—The graphs of M for h 2 = 0.3. Depicted for various values of the phenotypic deviation (in standard deviations), y 2 
varies from 0 to 0.7. 
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Mp^) = Pr(\C\ >K\P) = 1- 



'K 

F p,h*,Ag’ c ) dcd g■ 

-K 


(A17) 


Figure A1 depicts the graphs of Mas a function of y 2 for various absolute values of Pwhen /r > y 2 . Note that for a 
given A 2 , Prob(|C| > K\P) is always higher for a larger P , suggesting it is more productive to focus on classes of 
individuals with larger phenotypic deviations from the population mean. 

In the case where A 2 may be <y 2 we get a similar increase in probability for larger phenotypic values, as depicted in 
Figure A2. 

We can now utilize M to assess whether the variance components we have estimated through the familial covariances 
warrant further molecular investigation. Instead of basing our threshold on an arbitrary value of the epigenetic 
heritability, we decide on an epigenetic deviation threshold K, typically at one standard deviation (of P), and a 
probability threshold, T, and plug the epigenetic heritability, the genetic heritability, and a phenotypic deviation into 
M. If M= Prob(| C | > K\P) > T, then we gain confidence that the heritable epigenetic deviation, | C|, is large enough to 
warrant looking deeper into organisms with a phenotypic deviation larger than \P\. 
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